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The exact formulae for spectra of equilibrium diffusion in a fixed bistable 
piecewise linear potential and in a randomly flipping monostable potential 
are derived. Our results are valid for arbitrary intensity of driving white 
Gaussian noise and arbitrary parameters of potential profiles. We find: (i) 
an exponentially rapid narrowing of the spectrum with increasing height 
of the potential barrier, for fixed bistable potential; (ii) a nonlinear phe- 
nomenon, which manifests in the narrowing of the spectrum with increasing 
mean rate of flippings, and (iii) a nonmonotonic behaviour of the spectrum 
at zero frequency, as a function of the mean rate of switchings, for randomly 
switching potential. The last feature is a new characterization of resonant 
activation phenomenon. 

PACS numbers: 05.40-a, 05.10.Gg, 02.50. Ga 

1. Introduction 

Spectral densities of fluctuations provide an important tool to character- 
ize physical systems, because they can be measured directly in experiments. 
The investigations of spectra are useful to observe and analyze the inter- 
play between fluctuations, relaxation and nonlinearity which are inherent 
to real physical systems. This interplay ranks among the most challenging 
problems of modern nonlinear physics and forms the basis of well-known 
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nonlinear phenomena like stochastic resonance resonant activation |21, 
noise-enhanced stability [Hill], ratchet-effect [31111, etc. 

The exact formulae for spectra of fluctuations in nonlinear dynamical 
systems were first derived for thermal diffusion in fixed potentials. Caughey 
and Dienes [7] pioneered in applying analytical method based on Laplace 
transform of conditional probability density to the first-order system with 
y-shaped potential. Another approach has its origins in the expansion of 
probability density of transitions in terms of Fokker-Planck kinetic opera- 
tor eigenfunctions. This method was applied in [S] for obtaining correlation 
function of a bistable system with rectangular potential profile. We would 
also mention theoretical and numerical calculations reported in refs. 
concerning the spectra of underdamped double-well system driven by white 
Gaussian noise. In these papers the spectral peaks corresponding to stan- 
dard resonance and transitions between steady states have been revealed. 
Stationary spectra of fluctuations for monostable and bistable potential pro- 
files, by analog simulations of underdamped stochastic system driven by 
colored noise, have been experimentally obtained in ref. jJO]. The model of 
one-dimensional Brownian motion in singular potential like the potential of 
hydrogen atom was investigated in 11 . Authors detected some region of 
power spectrum with 1// frequency dependence. 

Despite a lot of work has been done to analyze spectra of fluctuations 
in the presence of one noise source, there is however lack of investigation 
on the so-called two-noise system spectra of fluctuations. A paradigmatic 
model is the overdamped Brownian motion in a randomly fluctuating po- 
tential. This model is being studied intensively in view of wide application 
in physics, chemistry and biology. However, an exact analytical results have 
been obtained only for escape rates, as mean first-passage times and life- 
times m and stationary probability distributions of Brownian motion 
jl3j . In this paper we report the exact calculations of diffusion spectrum for 
Brownian particle moving in fixed double-well potential and dichotomously 
switching linear potential. Our theoretical results, based on Markovian the- 
ory and on Laplace transform of conditional probability density, are valid 
for arbitrary intensity of driving white Gaussian noise and arbitrary param- 
eters of potential profiles. We find: (i) a narrowing of the spectrum with 
increasing height of the potential barrier for fixed potential; (ii) a narrowing 
of the spectrum with increasing mean rate of flippings, and (iii) a nonmono- 
tonic behaviour of the spectrum at zero frequency, as a function of the mean 
rate of switchings, for randomly switching potential. This last behaviour is 
a new characterization of resonant activation phenomenon . 
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2. Basic equations 

Let us consider an overdamped Brownian motion in a fixed potential 



U (x) described by Langevin equation 

dx dU (x) 
dt dx 



(1) 



where x{t) is the position of Brownian particle, ^(t) is a (5-correlated Gaus- 
sian noise with zero mean and intensity 2D. The Fokker-Planck equa- 
tion, or Smoluchowski equation [Til, for the conditional probability density 
W{x,t |xo,0) of Markovian random process x{t), corresponding to (0), is 



dW _ d 
dt dx 



dU (x) 
dx 



W 



+ D- 



dx"^ 



with initial condition 



W{x, |xo, 0) = 5 {x — xq) . 



(2) 



(3) 



Let us assume that a stationary regime exists, then the probabilistic flow 
equals zero at x ^ ±oo 



dW 

D— + U' ix) W 
ox 



0. 



^ = ±00 



(4) 



The correlation function of Brownian particle displacement x{t) in a sta- 
tionary state can be calculated as [T^ 



+ CX) 



K[r]= j xqWoo {xq) dxo j xW{x,T\xQ,Q)dx, (5) 

— oo — oo 

where Woo{x) is the stationary probability density (SPD) |14( llfij 



Tyoo(x) =C-e-^(")/^, 



C 



(6) 



To obtain the correlation function K [r] we need to solve the second-order 
partial differential equation @ using eigenfunction expansion |5| El • How- 
ever, as shown in j3Ej, the determination of SPD ^ together with the 
Laplace transform method are sufficient for calculating the spectral density. 
In fact from Wiener-Khinchin theorem we have 



+ 00 

S{uj) = — / K [t] cos (wr) dr = -Re [k [iu] ] 



(7) 
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where K [p] is Laplace transform of ^^[t]. By Laplace transforming 0, 
with initial condition 0, we obtain 

+ 1 Y]-pY = -S {x - xo) , (8) 

i.e. a second-order ordinary differential equation, where Y {x,xo,p) is the 
Laplace transform of conditional probability density 

+ 00 

Y{x,xo,p)= J e-P^W{x,t\xo,0)dt. (9) 



According to Eqs. ^ and © we solve © with boundary conditions 

= 0. (10) 



x=itoo 



By using Eqs. © and Q the Laplace transform K [p] of the correlation 
function is 

+ 00 +00 

K[p]= J xqWoo (xq) dxo J xY {x,xo,p)dx. (11) 

— oo — oo 

Then, after substitution of p = itu in Hll|l we can find the spectral density 
S{uj) from ((T)). Thus for calculating spectrum it will suffice to solve ordinary 
differential equation Q and make double integration, instead of solving 
partial differential equation ((21). To end we need the explicit expression of 
the internal integral in (jllj) . By multiplying both parts of Q on x and 
integrating it over the total area, using boundary conditions (|1U|) , we obtain 



1 
P 



xY {x,xo,p)dx = — [Y {oo,xo,p) - Y {-oo,xo,p)] 

p p 

— oo 



U' {x)Y {x,xo,p)dx. (12) 



-oo 



3. Fixed bistable potential 

Let us calculate the spectral density for symmetric double- well piecewise 
linear potential (see Fig. 1) 
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xo /- 



Fig. 1. Double-well piecewise linear potential. 



Substituting l(T^ in (ISjl. lfTU]) we obtain the following equation for the 
Laplace transform Y{x,xo,p) of conditional probability density 



DY" - ^ [sgn (x) ¥]' -pY = -6{x- xq) 
with the conditions at reflecting boundaries x = itL 



DY' - -sgn (x) Y 



0, 



x=±L 



(14) 



(15) 



where sgn(a;) is the sign function. Because of normalization condition for 
Y{x,xo,p) 

L ^ 

Y (x, xo,p) dx = -, 



equation ((T^ gives 
G(xo,p) 



^ + ^--[Y{L,xo,p)-Y i-L, xo,p)] 



p p^L p 



IE f 

— / Y {x,xo,p)dx. 



(16) 



-L 



To derive the function Y{x,XQ,p) we consider first xq > and solve 
homogeneous equation ((TH) in regions — L <x<0, 0<x<xo,xo<x<L 
separately. Then we apply the continuity conditions at the points x = 
and X = Xo 

D [Y' (+0, xo,p) - Y' (-0, xo,p)] -j[Y (+0, xo,p) + Y (-0, xo,p)] = 0, 
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D [V (xo + 0, xo,p) - Y' (xo - 0, xo,p)] = -1, 

y(+o,xo,p) =y(-o,xo,p), 

y(xo + 0,xo,p) =y(xo-0,xo,p). (17) 

Solving H14() in above-mentioned regions and taking into account the bound- 
ary conditions H15() we arrive at 



y (x,xo,p) 

where Ai^2 = 
gives 



ci 



C2 • e 



^i'' + C3 • e^2^, 



C4 



-L < X < 0, 
< X < xo, (18) 
< X < L, 



E±^/WT^pDL?) /{2DL). Substitution of (HHI) in JTH) 

(Cl - Ci) 



G{xo,p) 
2E 



xo E D 

— H 1 

p p^L p 



Ai 



+ 



pL\ 



1 ^ 



-ML 



(19) 



Calculating unknown constants ci and C4 from the continuity conditions 
((T7jl and substituting theirs in H19() we have 



G(xo,p) = h -^T H 2T" 

p p^L p'^L 



Aie" 



-AiL 



Ase" 



-A2L 



+ 



AiL 



p2L A2e-^2L_Aie-^i^' 

To obtain the function G{xq,p) in the region xq < we use symmetry 
considerations. Because of the symmetry of the potential U{x), the SPD (jEJ 
is an even function of x, W{x, t |xo, 0) = 1^(— x, 1 1— xq, 0) and Y (x, xo,p) = 
Y {—X, —xq,p). So G{xo,p) is an odd function of variable xq: G{xo,p) = 
—G{—xo,p), and from (|^. ((TT|) we obtain 

L 

/" „\ _/3a;n/Z/j„ (21) 



K[p] 



(e/^ -l)L 



xoG{xo,p) e^"«/^dxo, 



where (3 = E /D is the dimensionless height of potential barrier. Substitu- 
tion of (|2Uj) in (|21|) and subsequent integration gives the following result for 
Laplace transform of correlation function in stationary state 



k\p] 



(x2) 



D 

P 



,2 + p2 (1 



(3D 



+ 4/3 



e-l^) (aie"2 - ase^i) 
sinh2(a2/2) sinh^(ai/2)' 



02 



(22) 
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where ai^2 = i + 4pL^ /D^ /2. To obtain the spectral density of 
coordinate fluctuations of Brownian particle moving in a double- well poten- 
tial H13|) it remains to put in H22|) p = iuj and find its real part. However, 
we will not report here the exact formula for the spectrum because of its 
complicated expression. We give here the spectrum for particular case of a 
rectangular potential well, i.e. in the absence of a barrier (/3 = 0). We have 
(x^) = /3 and from we get 



so after substitution oi p = iuo in ()23() we arrive finally at {uo > 0) 



(23) 



Sr.) = ^ I sinhLyaJTD +sinLv^2^ 



4. Discussions 

Spectral densities of Brownian diffusion, obtained from H22|) . for differ- 
ent values of the potential barrier height are plotted in Fig. 2. As shown 




Fig. 2. Spectral density of Brownian particle displacement for different values of 
dimensionless height ([3 — E/D) of potential harrier: curve 1-13 — 2, curve 2 - 
/3 = 3, curve 3-/3=4, curve 4-/3 = 5. Parameters are L = 1, D — I. 

in Fig. 2, the spectral density has a maximum at zero frequency, which is a 
general property of Markovian random processes. We see also that the spec- 
trum S{uj) narrows very rapidly with increasing height of potential barrier, 
and its value at zero frequency increases fast. This nonlinear phenomenon is 
due to very rare transitions between steady states when the barrier is high 
with respect to the noise intensity [T^- Brownian particles therefore move 
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within a potential well for most of the time, and their displacements vary 
very slowly. As a result, the width of spectral density decreases. 

To verify this hypothesis we compare the behaviours of the spectral 
width and the mean rate of transitions between steady states as a function 
of potential barrier height. First we find the value of spectral density at 
zero frequency 5(0) = A' [0] /vr. Let us expand the function in power 
series on small parameter 02- Then we express this parameter in terms of 
small parameter p 

pL'^ p^L^ 
a2 ^ TTTT + 



f5D 

and after calculation of the limit p ^ 0, we get 

_ (/?-l)%2/3_(^3_3^2^4^_4).^/^_5 
^^^^-^ /34 (e/3-1) • ^^^^ 

The value ^(O) increases therefore as an exponential law ^(O) ~ e^ //?^ , 
with increasing height of potential barrier (3 and takes the finite value 
2L^ /{IbiiD) for /? = 0. This value corresponds to a diffusion in rectan- 
gular potential well (see ()24() 'l. The width of the spectral density with a 
maximum at zero frequency can be defined as JS] 

/ 2\ 

n = / 5(a;)da;/5(0) (26) 

The variance of Brownian particle position in stationary state from Eqs. © 
and ^ is 



{P^ - 2/? + 2) • e^ - 2 



/32 (e/3 - 1) ^^'^^ 

and increases monotonically from the value L^/S, which takes for /? — > 0, 
to the value L^, which takes for /3 — > 00, due to the finite area of diffusion. 
Substituting Eqs. and (HTJ in Eq. we obtain 



^ ttD {13^ -2(3 + 2) -e^ -2 



2L2 _ 1)2 e2/3 _ (^3 _ 3^2 + 4^ _ 4) . g/3 _ 5- 

By introducing correlation time similar to H26|l 

+00 

K [t] cLt /K [0] 



(28) 
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we find from Eqs. and 

_ ttS (0) _ 

~ 1^ ~ 2n' 

and equation 1)281) gives the exact correlation time for bistable potential, 
recently obtained in The spectral width decreases monotonically with 
increasing height of potential barrier from the value SvrZ) / (4L^) , taken for 
/3 — > 0, to zero, taken for /? — > oo. 

The mean rate of transitions, from one stable state to the other, can be 
determined through the mean first passage time (MFPT) to reach the top 
of barrier (x = 0) from the bottom of well (x = L), by solving the following 
differential equation |^ (x > 0) 

Dt"{x) + '^t' (x) = -1 

with boundary conditions: r' (L) = 0, r (0) = 0. After simple calculations 
we get for t{L) 

L2 e^-l-/3 

and for mean rate of transitions between two steady states 

n-^. (29) 
^^-L2 e/3-1-/?- ^^^^ 




Fig. 3. The spectral width and the mean rate of transitions between two steady states 
vs the dimensionless height of potential barrier. Parameters are L — 1, D = 0.5. 

In Fig. 3 we report the behaviours of 11 and as a functions of dimen- 
sionless height of potential barrier /3. The curves expressed by Eqs. (^5]) and 
H29|l practically coincide at large values of /3. Thus, the mean rate of tran- 
sitions is approximately the spectral width of Brownian particle coordinate 
fluctuations in a stationary state. 
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5. Randomly switching monostable potential 

Let us consider now two-noise nonlinear system, namely, one-dimensional 
overdamped Brownian motion in a fluctuating potential described by the 
following Langevin equation 

| = -?2iM + j(,), (30) 

where ^(t) is white Gaussian noise with zero mean and intensity 2D, <I>(x, t) = 
U{x) + ar]{t)x, U{x) is the same potential but without barrier (E = 0), 
and r]{t) is Markovian dichotomous noise switching with mean rate be- 
tween the values ±1. In other words, we analyze Brownian difl'usion in 
monostable potential with two randomly switching stable states near re- 
flecting boundaries at x = ±L (see Fig. 4). 



0(x,f) 




Fig. 4. Randomly switching monostable potential. 

Let us rewrite for our case the closed set of differential equations for 
probability density W{x,t), recently obtained in jJHl, in the diffusion inter- 
val {—L, L) 

dW _ dQ ^d^W 
dt dx dx"^ 

— = -2.Q + a— (31) 

with the following conditions at reflecting boundaries 

(W + aQ)^^^^ = 0, [DQ' + aW)^^^^ = {), (32) 

where Q{x, t) = {7]{t)5{x — x{t))) = W {x, t) {r]{t) \x (t) = x) is an auxiliary 
function We use the same method as for fixed potential. 
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According to Eqs. and (j32I) and initial conditions for the functions 
W{x,t), Q{x,t): W{x,0) = 6{x — xq), Q{x,0) = 0, we solve the following 
system of differential equations in the interval {—L, L) 

DY" + aZ' -pY = -6{x - xq), 

DZ" + aY' -{p + 2v)Z = {) (33) 

with boundary conditions 

{DY' + aZ)^^^^ = 0, {DZ' + aY)^^^^ = 0. (34) 

Here Y {x,XQ,p) and Z{x,xo,p) are the Laplace transforms of conditional 
probability density and of auxiliary function respectively. By putting E = 
in H16() we get 

G{xo,p) = ^ - - [Y{L,xo,p)-Y{-L,xo,p)] . (35) 
p p 

Now we solve the homogeneous set of linear differential equations (|33jl in 
two regions: —L < x < xq and xq < x < L. Then we find eight unknown 
constants from the boundary conditions ()34() and continuity conditions at 
the point x = xq 

^ |x=xo — — ^ |x=zo+0 ! ^ |x=2:o— — ^ |a::=a;o+0 ~l~ 1 / ; 

Z \x=xo—0 = Z \x=xo+0 , Z |a;=a;o-0 = Z |x=a;o+0 ■ 

After some algebra we obtain from ()35|) 

^(^0,^) = ^ (36) 
p 

{Dpi — p) sinh piL sinh p2Xo — {Dp2 — p) sinh P2L sinh PiXq 
p [p2 {Dpi — p) sinh piL cosh P2L — pi {Dp^ — p) sinh P2L cosh piL\ ' 

where 



Pl,2 



^t + :d^Vt + :d^' ^ = iD^^D- ^^^^ 



To find the Laplace transform 1)11(1 of correlation function in stationary 
regime we use the expression of SPD for our system, derived in ref. [TSj . 

1 1 + /i cosh 7x /cosh 7L 
^-("^ = 2L- l + ^tanh7L/(7L) ' ^^'^ 
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where /i = / {2vD) . After substitution of Eqs. (|H6j) and (|H8|) into Eq. 
((TT|l and integration we get 



where 



K[p] 



1 



p p[l + fitanhjL /{jL)] 
{Dpi -p) R (p2) tanh piL - {Dpi - p) R (pi) tanh P2L 
P2 {Dpi — p) tanh piL — pi {Dp^ — p) tanh P2L ' 




tanh 



(40) 



To obtain the exact formula for the spectral density of Brownian particle 
position it remains to put in equation (|39jl p = iuj and find the real part 
of expression. In the absence of flippings (a = 0, /i = 0) we find from Eq. 
(|H7j) : pi = yj^'^ + p /D , p2 = \fpjD and obtain the result for rectangular 
potential well of equation 



6. New characterization of resonant activation 

The evolution of spectrum shape with varying switchings mean rate 
V is shown in Fig. 5. The spectral density of Brownian diffusion in this 
non-Markovian case has also a maximum at zero frequency. For very large 




Fig. 5. Spectral density S (uj) for different values of switchings mean rate v: curve 
1 - v — 0.01, curve 2 - v = 3, curve 3 - v — 30. The curve 4 corresponds to a free 
diffusion in rectangular potential well {a — 0). The parameter set is: L — 1, a — 3, 
D = 0.5. 



spagnolo printed on February 2, 2008 



13 



values of z^, the spectral density approximates to the curve corresponding to 
a free diffusion in rectangular potential well. The main feature of Fig. 5 is 
that the spectrum at zero frequency S (0) shows nonmonotonic behaviour 
with increasing switchings mean rate v. Namely, S (0) initially decreases, 
reaches a minimum and then increases reaching asymptotically the value 
/ (ISttD) , obtained for rectangular potential well. 

Let us find the analytical expression of 5 (0) = K [0] /vr . Using the 
approximate expressions for parameters pi, p2 at small p (see (jHZI)) 

P ( 0^ \ \/2vT) 



and formula for the variance ^H] 



„2\ _ 7^^^ + 3^ [(2 + 7^^^) tanh7L - 27L] 



/ 373L[l + /itanh7L/(7L)] ^ ^ 

we get from Eq. 

S(0) = — -{IQuD-f^L"^ 

+50^ [60 + 27 fi + i'y'^L'^{3fi - 1) + {i'j'^L'^ - 27 fi - 12)-fLcothjL 

-(48 + (^ + 4) - 47^L4)^^^]}. (42) 

7L 

The typical z^-dependence of spectral density at zero frequency is plotted 
in Fig. 6. We see a clear minimum at z^ ~ 3. To explain this minimum let 
us consider the resonant activation phenomenon for this system. From the 




Fig. 6. Nonmonotonic behaviour of the spectral density at zero frequency as a func- 
tion of mean rate of flippings, for the same parameters L, a, D of Fig. 5. 
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closed set of differential equations for MFPTs T+ (x) and T_ (x) \1H\ 

DT'l -aT^ + u (T_ - T+) = -1, 
DT'l -aT'_ + u (T+ - r_) = -1 



(43) 



we calculate 7+ (x) and r_ (x), i.e. the MFPTs for positive (0) = +1 
and negative (0) = —1 initial value of the dichotomous noise, with start- 
ing position of Brownian particles at the point x, respectively. If we place 
the absorbing boundary at the point x = L we solve equations (|43j) with 
the following boundary conditions: Tj_ {—L) = 0, T-j. (L) = 0. The arith- 
metic average of MFPTs T (x) = [r+ (x) + T_ (x)] /2 for initial position of 
Brownian particles at the point x = —L is 



4zaL2 ^ 



cosh 27L — 1 



(sinh 27L - 27L)' 
cosh 27L -|- n 



(44) 



The behaviour of T[—L) as a function of switchings mean rate has a 
minimum, as shown in Fig. 7. This effect was called in literature reso- 



1 



T{-L) 



MFPT 



2 4 6 8 10 12 



Fig. 7. Resonant activation phenomenon for MFPT T {—L). The parameters L, a 
and D are the same as in Fig. 6. 

nant activation: the average residence time as a function of the barrier 
fluctuation rate v has a minimum at intermediate rates between very slow 
and very fast fluctuations |2]. In this range of rate the crossing event 
is strongly correlated with the potential fluctuations and Brownian parti- 
cles overcome randomly switching barrier in a minimal time. As a result, 
Brownian particle position changes rapidly and very slow components of 
the random process x {t) are present in minor amounts: the spectral density 
at zero frequency takes a minimum. Thus, the nonmonotonic behaviour 
of the spectral density at zero frequency S (0) can be interpreted as a new 
characterization of resonant activation phenomenon. 
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Finally we report in Fig. 8 the behaviour of spectral width H as a func- 
tion of flippings mean rate f. We find a new nonlinear phenomenon: the 




10 20 30 40 V 



Fig. 8. Spectral width vs flippings mean rate for the same parameters L, a and D 
of Fig. 5. 

spectral width decreases with increasing mean rate of switchings, contrary 
to the linear behaviour. As switchings mean rate increases, the slope of 
the potential profile of Fig. 4 becomes less and less important. As a re- 
sult, the diffusion time of Brownian particle between the reflecting bound- 
aries X = ibL increases and is determined by a free diffusion at very fast 
flippings. The random process x (t) therefore becomes more slow and the 
spectral width 11 decreases. 

7. Conclusions 

The exact formula for the spectral density of diffusion in double-well 
potential for arbitrary noise intensity and arbitrary parameters of poten- 
tial profile was obtained. We found very rapid narrowing of the spectrum 
with increasing height of a potential barrier between steady states. We 
also derived the exact result for spectral density of fluctuations in two-noise 
nonlinear system, namely, for overdamped Brownian diffusion in randomly 
flipping potential. We found a new characterization of resonant activation 
phenomenon in the behaviour of spectral density at zero frequency and new 
nonlinear effect associated with narrowing of the spectrum of Brownian 
particle position with increasing mean rate of switchings. Our analytical 
method enable us to investigate more difficult problems as those with more 
complex potential profiles. 

This work has been supported by INTAS Grant 2001-0450, MIUR, 
INFM, by Russian Foundation for Basic Research (project 02-02-17517), 
by Federal Program "Scientific Schools of Russia" (project 1729.2003.2), 
and by Scientific Program "Universities of Russia" (project 01.01.020). 
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